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THE SPATIAL ENERGY SPECTRUM OF MAGNETIC FIELDS IN OUR GALAXY 
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ABSTRACT 

Interstellar magnetic fields exist over a broad range of spatial scales, extending from the large Galactic scales 
(~ 10 kpc) down to the very small dissipative scales (<C 1 pc). In this paper, we use a set of 490 pulsars distributed 
over roughly one third of the Galactic disk out to a radius R ~ 10 kpc (assuming Rq = 8.5 kpc) and combine their 
observed rotation and dispersion measures with their estimated distances to derive the spatial energy spectrum of 
the Galactic interstellar magnetic field over the scale range 0.5- 15 kpc. We obtain a nearly flat spectrum, with 
a ID power-law index a = -0.37 ± 0.10 for E^{k) = Ck a and an rms field strength of approximately 6 /uG over 
the relevant scales. Our study complements the derivation of the magnetic energy spectrum over the scale range 
0.03- 100 pc by Minter & Spangler (1996), showing that the magnetic spectrum becomes flatter at larger scales. 
This observational result is discussed in the framework of current theoretical and numerical models. 
Subject headings: ISM: magnetic fields — pulsars: general — turbulence 



1. INTRODUCTION 

On the whole, the physical properties of the interstellar gas 
are better established and understood than those of the inter- 
stellar magnetic field. The spatial power spectrum of interstel- 
lar HI has been measured using observations of the 21 -cm line 
in several small areas across the sky. Emission measurements 
suggest a smooth power-law behavior with a 3D spectral index 
~ -3 for spatial scales > 5 pc (e.g., Crovisier & Dickey 1983; 
Green 1993; Dickey et al. 2001), while 21-cm absorption mea- 
surements yield a 3D spectral index ~ -2.75 in the scale range 
0.01-3 pc (Deshpande et al. 2000). 

The spatial structure of interstellar free electrons has been 
probed using observations of Galactic pulsars and extragalac- 
tic compact sources. Relying on interstellar scintillation data, 
Armstrong et al. (1995) showed that the spatial power spectrum 
of the interstellar free-electron density in the nearby interstel- 
lar medium (ISM) could be approximated by a single power 
law with a 3D spectral index ~ -3.7 (very close to the Kol- 
\ mogorov value of -1 1 /3) for scales ranging from ~ 10 10 cm to 



10 15 cm. By combining their scintillation data with rotation 
measure (RM) fluctuation measurements and with gradients in 
the average electron density, they were able to extend the range 
of their observed nearly-Kolmogorov spectrum up to a scale 
~ 10 20 cm(-30pc). 

Studies of pulsar and extragalactic radio-source RMs have 
yielded a wealth of observational information on the spatial 
structure of the Galactic interstellar magnetic field. However, 
most studies to date have focused on the large-scale field struc- 
ture (Simard-Normandin & Kronberg 1980; Han & Qiao 1994; 
Rand & Lyne 1994; Indrani & Deshpande 1998; Han et al. 
2002) and only a few have provided estimates for the rms am- 
plitude and the characteristic scale length of the turbulent mag- 
netic field. For instance, by adopting a single-cell-size model 
for the turbulent field and analyzing the residuals with respect 
to their best-fit model for the large-scale field, Rand & Kulkarni 
(1989) obtained a turbulent field strength of ~ 5 ^G and a cell 
size of ~ 55 pc. Ohno & Shibata (1993), who did not resort to 
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large-scale field models, organized their pulsar data by pairs of 
pulsars seen in nearly the same direction and interpreted them 
in the framework of a single-cell-size model for the random 
field. With an assumed cell size in the range 10-100 pc, they 
obtained a random field of amplitude 4-6 /iG 

As Rand & Kulkarni (1989) themselves acknowledged, the 
turbulent magnetic field cannot be satisfactorily described by 
a single scale length. Simonetti et al. (1984) and Simonetti & 
Cordes (1986) were the first authors to look into the question 
of how the turbulent magnetic energy in the ISM is actually 
distributed in the scale range ~ 0.01 - 100 pc. They used the 
structure function of RMs of extragalactic sources to investi- 
gate variations in the composite quantity n e B, where n e is the 
free-electron density and B is the magnetic field. Their study 
suggests that the 3D spatial power spectrum of fluctuations in 
n e B at scales > 0.01 pc can be described by a power law of 
index ~ -3.1 with an outer scale estimated at < 90 pc. Later, 
Minter & Spangler (1996) performed a combined analysis of 
the structure functions of RM and emission measure (EM) to 
separate the fluctuations in electron density and magnetic field 
in the scale range ~ 0.03 - 100 pc. They measured the RMs of 
38 extragalactic sources located in a small area of the sky near 
Galactic coordinates / = 144°, b = 21°, where the EMs were 
deduced from measured Ha intensities. For scale sizes up to 
~ 4 pc, the structure functions of RM and EM were consistent 
with a 3D-turbulence model in which both electron density and 
magnetic field fluctuations have Kolmogorov spectra. Between 
~ 4 pc and ~ 80 pc, the structure functions were consistent with 
2D turbulence. Minter & Spangler (1996) also obtained an rms 
amplitude of ~ 1 /iG for the turbulent magnetic field on scales 
up to 4 pc. 

It is clear that the magnetic energy spectrum does not abruptly 
stop at the estimated outer scale of the turbulence, ~ 80 pc. The 
spectral energy of the magnetic field will extend continuously 
up to the largest Galactic scales, as a result of an inverse cas- 
cade of magnetic helicity (e.g., Pouquet et al. 1976) and/or due 
to other physical processes involving compression and shear- 
ing at large scales (e.g., Parker 1971; Kulsrud 1986). turbulent 
field and the more structured field at larger scales describe the 
present-day properties of the Galactic magnetic field and pro- 
vide the necessary observational reference for theoretical mod- 
els. Knowledge of the complete magnetic energy spectrum can 
offer a solid observational test for dynamo and other theories 
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FIG. 1 . — Distribution projected on the Galactic plane of the 490 pulsars with well-determined RMs viewed from the North Galactic pole. Pulsars with a positive 
RM are denoted with crosses, those with a negative RM are denoted with circles, and the symbol area is proportional to the RM. Superimposed onto the pulsar 
distribution is an approximate outline of the four known spiral arms. 



for the origin of Galactic magnetic fields. 

The purpose of the present paper is to determine the por- 
tion of the magnetic spectrum that lies above the strictly tur- 
bulent domain, thereby filling the observational gap between 
the large Galactic scales and the small turbulent scales. Pulsars 
are unique in providing a direct measure of the average line-of- 
sight component of the interstellar magnetic field, weighted by 
the local electron density, from the ratio of RM and dispersion 
measure (DM). Here we will analyse the available pulsar RM 
and DM data in combination with their estimated distances to 
deduce the interstellar magnetic energy spectrum E B across the 
scale range 0.5 - 15 kpc. In § 2, we present our pulsar data and 
explain how they are used to construct the interstellar magnetic 
energy spectrum, in § 3, we review the assumptions underlying 
our derivation, and in § 4, we discuss the significance of our 
results and compare them with theoretical predictions. 

2. THE MAGNETIC ENERGY SPECTRUM FROM PULSAR DATA 
2. 1 . Derivation of (Eg )(S) 

At the present time, we have 521 pulsars with measured RMs 
at our disposal: 330 of them are published measurements mostly 
from Hamilton & Lyne (1987); Rand & Lyne (1994); Qiao et al. 
(1995) and Han et al. (1999), while another 202 RMs come 
from unpublished observations with the Parkes telescope (see 
Han et al. 2002). We discard all pulsars with a RM uncertainty 
larger than 30 rad m~ 2 , leaving us with a set of 490 pulsars. 
The distances to the pulsars are determined from their DM us- 
ing the NE2001 model for the Galactic electron density distri- 
bution (Cordes & Lazio 2003). The spatial distribution of our 
490 pulsars is displayed in Figure 1, showing that they cover 
roughly one third of the thick ISM disk out to a Galactocentric 
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FIG. 2. — Geometry of paired pulsars Pi and P 2 . indicates the position of 
the Sun and P[ represents the projection of Pi onto the line of sight to P2. 



radius R ~ 10 kpc. 

The average line-of-sight component of the interstellar mag- 
netic field between a pulsar (P) and the Sun (0), weighted by 
the local free-electron density, is given by the ratio of the pul- 
sar's RM to its DM: 

joVfiil dl 

"II \r. ~ - 



- (1.232 ^G) — 
J D n e dl DM 



(1) 



Here, n e is the free-electron density (in cm" 3 ), B\\ is the line-of- 
sight component of the magnetic field (in ^G), dl is the length 
element along the line of sight to the pulsar (in pc), D is the 
distance to the pulsar (in pc), and DM and RM are expressed 
in pc cm -3 and in rad m~ 2 , respectively. Likewise, the electron- 
weighted average value of B\\ between two pulsars, Pi and P2, 
located on the same line of sight is given by 

fin| Pi =(1.232 t**) DM2 _ DMi - (2) 

In order to extract as much information as possible on the 
spatial variations of the interstellar magnetic field, we are go- 
ing to use both equation (1) applied to individual pulsars and 
equation (2) applied to pairs of pulsars lying in almost the same 
direction. Here, two pulsars Pi and P2 are considered to lie in 
almost the same direction if their angular separation, 6, satisfies 

d|pj=Di sin6»<0.1 min(D b D 2 -Di) , (3) 
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FIG. 3. — Measured values of \B« | as a function of 5, where 5 is either 
the pulsar distance or the separation of pulsar pairs. The 490 filled circles 
correspond to measurements to individual pulsars while the 1200 open circles 
correspond to measurements between paired pulsars. 



supposing that Pi is the closer pulsar (see Figure 2). As we 
will see from equation (13) and Figure 5 below, B|||p* is only 
sensitive to magnetic fluctuation cells with sizes > 0.15 (where 
S = D\ orD 2 -D\). Equation (3) ensures that line segments ©Pi 
and ©Pj pass through the same magnetic fluctuation cells with 
sizes > 0.15. The maximum angular separation between two 
paired pulsars is 5.7°. Note that our pair-selection approach is 
similar to that originally proposed by Ohno & Shibata (1993) 
for estimating the rms amplitude of the random field. 

Altogether, we have 490 measurements of to individual 
pulsars and 1200 measurements between paired pulsars. For 
the following, we denote by 5 the distance over which B» is 
measured, i.e., S = D for the first type of measurements and 
5 = L>2 —D\ for the second type. Measured values of | are 
shown as a function of 5 in Figure 3. A few pulsar pairs had 
similar DMs but very different RMs, resulting in anomalously 
high values of By. In fact, the separation of these pulsars is 
almost certainly greater than indicated by the DMs, so these 
values were capped at 15 fiG. These and the few values between 
10 and 15 /iG are plotted at the top of the figure. 

The magnetic energy associated with Bu for each pulsar or 

2 

pulsar pair, Eg =B|| /(87r), is plotted as a function of ks = 1 /5 
in Figure 4 on logarithmic scales. There is a wide scatter in the 
points, largely because, for any pulsar located at a given dis- 
tance 5 and for any pair of pulsars separated by a given 5, \Bn | 
and Eg can vary from zero to some maximum value depend- 
ing on the angle between the line of sight and the direction of 
the local magnetic field as well as on the position of the pulsar 
or pulsar pair with respect to the maxima and minima of the 
magnetic fluctuations. Never-fhe-less, by averaging Eg over 
successive intervals along the fc s -axis, we can obtain estimates 
for the mean value of the magnetic energy associated with Bp , 

^ , (4) 

as a function of k s . It is important to realize that the above aver- 
aging procedure (denoted by angle brackets) implies averages 
over many different Galactic locations and many different lines 
of sight. Despite the wide scatter in Eg , the average values ob- 
tained for (Eg } (indicated with thick crosses in Figure 4) show 




a clear trend with ks, at least for k s < 2 kpc 1 . A least-squares 
fit of a power law to our estimated (Eg } gives 



(Eb u )(S) = Co 







(5) 



with [3 = 0.66 ±0.10, Co = (5.1 ±0.7) x 10" 13 erg cm" 3 , k s = 1 /5, 
and the normalising factor ko = 1 kpc -1 . 

Individual pulsar distances are uncertain by ~ 20% or about 
0.5 kpc at distances of a few kpc and there are only few pul- 
sars within 0.5 kpc from the Sun. Therefore, our estimates for 
the parallel field strength and energy are unreliable at scales 
< 0.5 kpc. In addition, there are only a few pulsars with dis- 
tances greater than 15 kpc, so the range of validity of equa- 
tion (5) is approximately 0.5 kpc < 5 < 15 kpc. 

2.2. Relationship between (Eg )(5) and E B (k) 

The next important step is to relate our measured (Eg )(5) 
to the magnetic-energy spectral density, E B (k). We proceed 
on the assumption (discussed in § 3) that the spatial fluctua- 
tions of the magnetic field, B, are statistically homogeneous 
and isotropic. We start by writing B(r) in terms of its three- 
dimensional Fourier transform, B(k): 
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B(r)= B(k)e 



„27rik-r 



dk 



(6) 



where r is the position vector and k is the wave vector defined 
such that the relationship between wavenumber, k, and wave- 
length, A, is k = 1 /A. The line-of-sight component of B at a 
distance / can then be written as 
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Bu(k)e 27Tik tt' dk' 



(7) 



Its average value between two points located on the same line 
of sight at distances D\ and D2, respectively, reads 

r D 2 



— 1 



B||(/) dl 



D2-D1 JDl 

or, after substitution of equation (7) and integration over I, 
Bii(k) sinc(7rfe N 5) e 27 "'*" Dm dk , 



B 



(8) 



(9) 



with S = D2—D1, D m = (Di +Z>2)/2, and sinc(x) = (sinx)/x. 

The expectation value of the magnetic energy associated with 
Bp may be defined as 



8tt 
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8^ 



(B||(k) B|(k )} sinc(7r/fc|| S) sinewy S) 
e 2ni(k rV D " dkdk ■ (10) 



In the last line of the above equation, the expectation value is 
applied to the factor B||(k) Bjj(k ), which is the only stochastic 
part of the expression. Now, for homogeneous and isotropic 
turbulence, we have 
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FIG. 4. — Measured values of the magnetic energy associated with By, Eg =B|| /(87r), versus wavenumber, ks = 1/5, where S is either the pulsar distance or the 

separation of pulsar pairs. The thick crosses give average values (£ 5 } and their uncertainties over the 7 successive intervals delimited by the vertical dotted lines 

and corresponding, from right to left, to S(kpc) = [0.125,0.25], [0.25,0.5], [0.5,1.0], [1.0,2.0], [2.0,4.0], [4.0,8.0], [8.0,16.0]. The oblique thick line is a power-law 
fit to the {Eg ) with k s < 2 kpc -1 (see eq. [5]). The number of data points in each interval is indicated at the bottom. 



where £fi(k) is the 3D magnetic-energy spectral density and 
E B (k) is its ID counterpart. In view of equation (11), equa- 
tion (10) reduces to 

(E g|| )(S) = I J ^ sinc 2 (7rfe|| S) dk • (12) 

The integration over the direction of k is conveniently carried 
out in a cylindrical coordinate system with polar axis along the 
line of sight, such that hi = k cost? and dk = 2nk 2 dk d(cos9). 
The final result can be cast in the form 

(\XS)=\J E B (k)y(kS)dk, (13) 

with 

1 r ks /sinx\ 2 

y(kS)=— / ) dx. (14) 

nkS J a \ x J 

Clearly, the 1/3 prefactor in equation (13) follows from our as- 
sumption of isotropy whereby one third of the magnetic energy 
goes into the line-of-sight component. 

The shape of the weight function, y(kS), is displayed for dif- 
ferent values of the pulsar separation S in Figure 5. For any 
given S, magnetic fluctuations with kS <C 1, i.e., with wavenum- 
ber k <C ks or wavelength \^>S, give their full weight to (Eg )(S) 
Fluctuations with k ~ ks or A ~ S give approximately half their 
weight, and fluctuations with k » k s or A -C S make a neg- 
ligible contribution. Hence the main contribution to (Eg )(5) 

comes from fluctuations with wavenumber k < ks- Physically, 
for long-wavelength fluctuations (A » S), B\\ is nearly constant 
between both pulsars, so that its average value B» is approxi- 
mately equal to its local value. For short-wavelength fluctua- 
tions (A <C S), B|| reverses many times between both pulsars, so 
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FIG. 5. — Weight function, y(kS) (defined by eq. [14]), as a function of 
wavenumber, k, for 5 selected values of the pulsar separation, S (corresponding 
to the averaged S in the first 5 intervals of Figure 4), in the 3D case. The dashed 
lines show the weight function in the 2D case (eq. [16]). 

that the contributions to B|| from line segments with opposite 
sign of Bu cancel each other out almost entirely. For fluctua- 
tions with A = S, there is complete cancellation if k is along the 
line of sight, and there is no cancellation if k is perpendicular 
to the line of sight; an average over the direction of k then leads 
to a net weight factor ~ 0.5. 

2.3. Determination ofE B (k) 

In practice, we tried various analytical functions for E B (k). 
For each of them, we computed the corresponding function 
(Eg }(S) by integrating the right-hand side of equation (13) nu- 
merically and then compared the computed (Eg ) to the mea- 
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sured } displayed in Figure 4. The lower limit of integra- 
tion as well as the step length of k were set to 0.02 kpc" 1 . The 
upper limit of integration is mainly determined by the shape of 
the weight function, y(kS), shown in Figure 5; here we adopted 
the value at which the function E B (k)y(kS) drops below 0.1% 
of its peak value. 

The trial function E B (k) that provides the best fit to the mea- 
sured (Eg ) is a slowly decreasing power law described by 



E B (k) = C 



k 



(15) 



with a = -0.37 ± 0. 10 and C = (6.8 ± 0.8) x 10" 13 erg cm" 3 kpc. 
As a reminder, k is the wavenumber defined as 1 / A and ko = 
1 kpc" 1 . The power-law indices of (Eg )(S) (see eq. [5]) and 
E B (k) (see eq. [15]) are approximately related through (3 = 1 +a. 
This relation would be exact if the weight function y(kS) were 
approximated by the step function u(l-kS), where u(x) is de- 
fined such that u(x) = for x < and u(x) = 1 for x > 0. Con- 
sequently, the uncertainty in a is comparable to the uncertainty 
in (3, i.e., ~ 0.10. Furthermore, since (Eg )(S) has been reli- 
ably determined over the range 0.5 kpc < S < 15 kpc (see Sec- 
tion 2.1), the validity range of equation (15) is approximately 
0.07 kpc" 1 <Jfc<2kpc _1 . 

The rms value of the fluctuating magnetic field with wavenum- 
bers in the above range, obtained by integrating equation (15), 
is 6.1 ±0.5 fiG. We note that this rms field strength is consis- 
tent with estimates by previous authors, e.g., via modelling dif- 
fuse 7-ray data and synchrotron emission (Strong et al. 2000), 
via energy equipartition (Berkhuijsen, referenced in Beck et al. 
1996), or based on pulsar data and synchrotron data by Heiles 
(1996). This leads us to believe that equation (1) using pulsar 
RM and DM data does not significantly underestimate or over- 
estimate the field strength in the ISM, at least on average over 
a large sample. 

3. VALIDITY OF THE UNDERLYING ASSUMPTIONS 

We have used the tools of homogeneous and isotropic turbu- 
lence theory to derive the spatial power spectrum of magnetic 
energy in our Galaxy. However, strictly speaking, these are not 
directly applicable to the problem at hand. Firstly, the spatial 
variations of the magnetic field at the scales probed here are 
not all stochastic; some are associated with large-scale coherent 
features such as supernova remnants or superbubbles and oth- 
ers are associated with the large-scale Galactic spiral structure. 
In this study, we did not try to distinguish the purely stochastic 
fluctuations from the coherent variations; all are included in our 
magnetic energy spectrum. This procedure is justified if the re- 
sults are taken to represent present-day Galactic averages rather 
than true statistical averages. 

A second, more fundamental, reason is that the magnetic en- 
ergy distribution is neither homogeneous across the sampled 
region (it includes spatial variations up to scales approaching 
the size of the region, which is itself a substantial portion of 
the Galaxy) nor isotropic (the effects of the disk geometry and 
the associated vertical stratification are significant at the rele- 
vant spatial scales). The most obvious manifestation of the lack 
of isotropy is the dominant horizontal component of the large- 
scale magnetic field. 

The failure of the assumption of homogeneity is inherent in 
our scale range. Minter & Spangler (1996) did not encounter 
the same problem because they studied only small-scale fluctu- 



ations, such that they were able to perform statistically mean- 
ingful spatial averages over a region which is small enough for 
turbulence to be assumed homogeneous across it. The down- 
side is that their spectrum applies only to their small region of 
observations and to small scales. In contrast, our derived spec- 
trum does not constitute a good statistical average. However, 
since our pulsar pairs are reasonably well distributed through- 
out the sampled volume, our spectrum may be regarded as a 
spatial average over the sampled volume - with some bias to- 
ward the vicinity of the Sun. Since the sampled volume corre- 
sponds roughly to a 120° wedge of the Galactic disk seen from 
the Galactic center, it is probably representative of the whole 
disk. In this case, our derived spectrum also provides a good 
estimate for the Galaxy-wide average spectrum. 

The assumption of isotropy is almost unavoidable and rou- 
tinely made in practice (e.g. Ohno & Shibata 1993; Minter & 
Spangler 1996), even though it fails both at large scales (for the 
reason outlined above) and at small scales (due to the presence 
of the large-scale magnetic field). Here, the use of the isotropy 
assumption and of the ensuing formulae would be justified if 
the sky coverage in line-of-sight direction were complete and 
uniform and if the magnetic spectrum were homogeneous (so 
as to give all directions the same weight). However, neither 
condition is fully satisfied. The resulting uncertainty can be es- 
timated once it is realized that the situation we are dealing with 
is intermediate between 3D isotropic and 2D isotropic (with the 
magnetic field at all locations and all sightlines being parallel 
to the Galactic plane). The 2D case can be treated similarly to 
the 3D case studied in § 2.2 and § 2.3: the expression found 
for (Eg^ )(S) is identical to equation (13) with the 1/3 prefactor 

replaced by 1/2 and the weight function given by 



y(kS) : 



1 



sinx 



dx . 



(16) 



n Jo ^(TrkS) 2 -x 2 
A numerical fit to Figure 4 then leads to a corresponding equa- 
tion (15) with a = -0.39 ±0.10 and C 2D = (8.6 ± 1.0) x 10" 13 
erg cm" 3 kpc. This shows that the derived values are not very 
sensitive to the isotropy assumption. The reason why a sim- 
ilar spectral index is obtained in the 2D and 3D cases is eas- 
ily understood. In § 2.3, we saw that the spectral index ob- 
tained with the step-function approximation to the weight func- 
tion (a = [3- 1 = -0.34) is close to that obtained with the correct 
weight function (a = -0.37). This indicates that the slope of the 
energy spectrum depends little on the exact shape of the weight 
function, but is governed by the relative positions of the cutoff 
wavenumbers. 

Another important assumption implicit in our derivation is 
the absence of correlation between the fluctuations in electron 
density and those in magnetic field strength, which allows us to 
approximate the space-averaged value of Bp between two points 
(eq. [8]) by a density-weighted average (eq. [1] or [2]). This as- 
sumption may not be well satisfied in the ISM. Both positive 
and negative correlations are possible and their net effect is dif- 
ficult to quantify (Beck et al. 2003). However, as mentioned 
in § 2.3, our derived rms field strength is consistent with other 
estimates. In addition, the fluctuation spectral index is affected 
only if the degree of correlation is scale-dependent. 

Finally, individual pulsar distances are uncertain by a large 
factor, typically ~ 20%, although the overall distance scale is 
more accurate. The impact on our derived magnetic spectrum 
is limited by the fact that the values of (Eg )(5) were obtained 
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from the average of many measurements of B» . 
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FIG. 6. — Composite magnetic energy spectrum in our Galaxy. The thick 
solid line is the large-scale spectrum derived in this paper. The thin solid 
and dashed/dotted lines give the Kolmogorov and 2D-turbulence spectra, re- 
spectively, inferred from the Minter & Spangler (1996) study. Their Kol- 
mogorov spectrum, after adjusting to our definition of k, can be written as 
E B (k) = C K (k/k r 5/3 with C K = 9.5 X 1CT 13 erg cnT 3 kpc and ko = 1 kpc" 1 . 
The 2D-turbulence spectrum is uncertain (see main text); it probably lies be- 
tween the dashed (Eb(&) <x k~ 2 / 3 ) and dotted (£bW oc At 5 / 3 ) lines, which are 
both scaled to match the Kolmogorov spectrum at £ 3D = (4 pc)~' . 



4. DISCUSSION 

We have derived the interstellar magnetic energy spectrum 
over the scale range 0.5- 15 kpc, based directly on measure- 
ments of pulsar RMs and DMs. The outer scale is a significant 
fraction of the size of the Galaxy and so contributions from the 
so-called large-scale (or regular) magnetic field are included. 
The spectral region covered by our derived spectrum has never 
been studied observationally before, yet it constitutes a crucial 
piece of the complete spectrum, whose knowledge is important 
for our understanding of the behavior and evolution of Galactic 
magnetic fields. 

In Figure 6, we plot our measured spectrum with that inferred 
from the Minter & Spangler (1996) study at smaller scales. This 
spectrum has a Kolmogorov slope (E B (k) oc fc" 5 / 3 ) up to a scale 
l 3D rv> 4 pc and an uncertain shape consistent with 2D turbulence 
between lj D and I2D ~ 80 pc. If, as suggested by Minter & Span- 
gler (1996), the 2D turbulence is in the form of thin sheets with 
thickness ~ /3D, we could assume that fluctuations with scales 
larger than l^u persist only if their projected wavelength normal 
to a sheet, A cos#, is less than the sheet's thickness lj D . The 
fraction of acceptable wave vectors at a given wavenumber k 
is then hv/\ = k/k^ D . If the acceptable wave vectors preserve 
a Kolmogorov spectral energy, the magnetic spectrum in the 
wavenumber range &2d - &3D is simply the Kolmogorov spec- 
trum multiplied by the above factor, i.e., a k~ 2 / 3 spectrum. On 
the other hand, if 2D turbulence means that the spectral energy 
is simply redistributed without loss into two dimensions, then 
the Kolmogorov spectrum oc A:" 5 / 3 extends out to kzo- The ac- 
tual situation probably lies between these two extremes. 

The two important features emerging from the composite 
spectrum (Figure 6) are (1) a change in slope from a = -5/3 at 



the small (< 4 pc) scales to a ~ -0.37 at the large (> 0.5 kpc) 
scales and (2) a possibly significant discontinuity between ~ 
80 pc and ~ 0.5 kpc. How do these two features fit with theo- 
retical predictions? 

According to our theoretical understanding of the ISM, tur- 
bulent energy is injected into the ISM by stellar winds and su- 
pernova explosions on scales ~ 10-100pc. Magnetic energy is 
then transferred, via a so-called direct cascade, toward smaller 
and smaller scales. In addition, because interstellar turbulence 
is helical, there is also an inverse cascade of magnetic helicity 
toward larger and larger scales, making it possible to amplify or 
maintain a large-scale magnetic field (Pouquet et al. 1976). Be- 
cause the quantity that is transferred as a whole in the domain 
of the inverse cascade is magnetic helicity (Hs(k) oc E^,(k)/k) 
rather than magnetic energy, the corresponding spectrum is flat- 
ter than in the domain of the direct cascade. Therefore, our 
finding a flatter spectrum at the larger scales is consistent with 
the inverse cascade theory. 

However, the value obtained for the spectral index (a ~ -0.37) 
differs from current theoretical estimates. In the theory of 3D 
MHD helical turbulence developed by Pouquet et al. (1976), a 
numerical solution of the evolution equation for the magnetic 
energy spectrum yields Es(k) oc k~ l at scales larger than the in- 
jection scale. As shown by these authors, the same spectral 
slope is obtained by applying a dimensional Kolmogorov-type 
argument to the inverse cascade of magnetic helicity. Direct 
numerical simulations do not yet convey an unambiguous or 
consistent picture. Balsara & Pouquet (1999) obtained numer- 
ical spectra consistent with the k~ l theoretical prediction. In 
contrast, in the simulations of Brandenburg (2001), a wave of 
spectral magnetic energy propagates toward smaller k under a 
k~ l envelope, but no k~ l spectrum actually develops. Instead, 
magnetic energy builds up at the largest scales and a secondary 
peak appears at the injection scale. Maron & Blackman (2002) 
showed that the actual shape of the magnetic energy spectrum 
depends in fact quite sensitively on the level of fractional ki- 
netic helicity with which the turbulence is forced. According to 
their simulations, high levels of kinetic helicity lead to double- 
peaked spectra similar to that obtained by Brandenburg (2001), 
whereas low levels of kinetic helicity give rise to a single peak 
at the small resistive scale. It is clear that one has to be ex- 
tremely cautious when trying to apply numerical spectra to the 
real Galactic magnetic field. Current simulations are still highly 
idealized (e.g., uniform background medium, incompressible or 
isothermal gas, periodic boundary conditions), and they neglect 
potentially important factors such as the true nature of turbu- 
lent forcing, the presence of large-scale shear, and the actual 
boundaries of the Galactic disk. Another difficulty is that the 
saturated spectra obtained numerically emerge after resistivity 
comes into play, whereas the Galaxy is probably not yet in the 
resistive regime (Blackman, private communication). 

The discontinuity appearing between the <100pc and >0.5kpc 
scales in the composite spectrum of Figure 6 can be partly ex- 
plained by the fact that the different parts of the spectrum apply 
to different Galactic regions. As discussed in § 3, our portion of 
the spectrum was obtained from a magnetic energy average over 
about one third of the inner Galactic disk, while that of Minter 
& Spangler (1996) was obtained from a small high-latitude re- 
gion in the outer Galaxy, where the fields are presumably much 
weaker than in the inner disk. However, if the 2D-turbulence 
portion of the spectrum is as flat as k~ 2 ^, the observed disconti- 
nuity is so severe that it suggests a genuine excess in the Galac- 
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tic magnetic energy spectrum at large scales. This excess com- 
pared to the amount expected from the inverse cascade alone 
could indicate that most of the energy input to the large-scale 
magnetic field occurs directly at large scales. This interpre- 
tation would be consistent with the standard Galactic dynamo 
theory, in which the large-scale magnetic field is amplified or 
maintained through the combined action of small-scale helical 
turbulence, presumably via the inverse cascade leading to the 
so-called alpha-effect, and large-scale shear, typically associ- 
ated with the Galactic differential rotation (e.g., Parker 1971; 
Vainshtein & Ruzmaikin 1971). It is widely believed that the 
large-scale shear acts much more efficiently than the helical tur- 
bulence (e.g., Parker 1971; Ruzmaikin et al. 1988; Brandenburg 
et al. 1992; Ferriere & Schmitt 2000; Blackman 2000), espe- 
cially when the spiral structure of the Galaxy and the stream- 
ing motions along spiral arms are taken into account (e.g., Ro- 
hde et al. 1999; Elstner et al. 2000). A more extreme point of 
view would be to consider that the large-scale magnetic field is 
not affected at all by the small-scale turbulent field and that all 
the energy contained in the large-scale field is directly injected 
at large scales. Various scenarios have been proposed in this 
spirit, such as protogalactic collapse and subsequent shearing 
by the Galactic differential rotation with rapid ambipolar dif- 
fusion (Kulsrud 1986) or protogalactic collapse and subsequent 
excitation of spiral arms and bars resulting in nonaxisymmetric, 
not-perfectly-azimuthal motions (Lesch & Chiba 1997). 
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